%%bash
ls /u/ki/swmclau2/des/SherlockPearceMCMC/*.hdf5 -ltr
Out[353]:
['/u/ki/swmclau2/des/SherlockPearceMCMC/VpeakSHAM_xi_gg_rmin_1.0_CorrAB.hdf5',
'/u/ki/swmclau2/des/SherlockPearceMCMC/VpeakSHAM_xi_gg_rmin_1.0_HOD.hdf5',
'/u/ki/swmclau2/des/SherlockPearceMCMC/VpeakSHAM_xi_gg_rmin_1.0_HSAB.hdf5']
Out[355]:
[u'emu_type',
u'training_file',
u'fixed_params',
u'emu_hps',
u'true_data_fname',
u'true_cov_fname',
u'sim',
u'obs',
u'cov',
u'nsteps',
u'param_names',
u'nwalkers',
u'nburn',
u'seed',
u'chain_fixed_params']
Out[361]:
array(['ombh2', 'omch2', 'w0', 'ns', 'ln10As', 'H0', 'Neff',
'mean_occupation_centrals_assembias_corr1',
'mean_occupation_satellites_assembias_corr1', 'logM1', 'logM0',
'sigma_logM', 'alpha'], dtype='|S42')
if model == 'HOD':
param_names = [r'$\Omega_b h^2$', r'$\Omega_c h^2$', r'$w_0$', r'$n_s$', r'$\ln(10A_s)$', r'$H_0$', r'$N_{eff}$',\
r'$\log(M_0)$', r'$\sigma_{\log M }$', r'$\log(M_1)$', r'$\alpha$']
elif model == 'HSAB':
param_names = [r'$\Omega_b h^2$', r'$\Omega_c h^2$', r'$w_0$', r'$n_s$', r'$\ln(10A_s)$', r'$H_0$', r'$N_{eff}$',\
'$\mathcal{A}_{sat}$',r'$\mathcal{A}_{cen}$', r'$\log(M_1)$',r'$\log(M_0)$', r'$\sigma_{\log M }$',r'$\alpha$' ]
else: #CorrAB
param_names = [r'$\Omega_b h^2$', r'$\Omega_c h^2$', r'$w_0$', r'$n_s$', r'$\ln(10A_s)$', r'$H_0$', r'$N_{eff}$',\
'$\rho_{sat}$',r'$\rho_{cen}$', r'$\log(M_1)$',r'$\log(M_0)$', r'$\sigma_{\log M }$',r'$\alpha$' ]
Out[366]:
<chainconsumer.chainconsumer.ChainConsumer at 0x7fa07714ad90>
Out[367]:
{'gal_type': 'SHAM',
'halo_property': 'halo_vpeak',
'min_ptcl': 100,
'nd': '5e-4',
'scale_factor': 1.001,
'scatter': 0.2,
'sim_hps': {'system': 'sherlock'},
'simname': 'ds_14_b_sub'}
cosmo_true_vals = [sim_info['cosmo_params'][key] for key in chain_param_names if key in sim_info['cosmo_params']]
hod_true_vals = [sim_info['hod_params'][key] for key in chain_param_names if key in sim_info['hod_params']]
/u/ki/swmclau2/.conda/envs/hodemulator/lib/python2.7/site-packages/matplotlib/figure.py:403: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
"matplotlib is currently using a non-GUI backend, "
fig = c.plotter.plot(figsize=(10,10), parameters = cosmo_param_names, truth = cosmo_true_vals)# parameters = [param_names[i] for i in (1, 4)])
#, truth = true_vals)
fig.show()
/u/ki/swmclau2/.conda/envs/hodemulator/lib/python2.7/site-packages/matplotlib/figure.py:403: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
"matplotlib is currently using a non-GUI backend, "
/u/ki/swmclau2/.conda/envs/hodemulator/lib/python2.7/site-packages/matplotlib/figure.py:403: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
"matplotlib is currently using a non-GUI backend, "
c.plotter.plot_walks(parameters = cosmo_param_names, truth = cosmo_true_vals)